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D . Abstract 

D , We report on preliminary results of searches for penguin mediated B decays based on 20.7 fb ^ 
' of data collected at the T(4S) peak with the BABAR detector at PEP-II. The following branching 
fractions have been measured: B{B+ (l)K+) = (7.7tJ-^ ±0.8) x 10-^ ^3(5° ^ cpK^) = {S.ltH ± 
V : 0.8) X 10-^ B{B+ (t)K*+) = (9.7tt^ ± 1.7) x 10-*^^ B{B^ (pK*^) = (8.7l|f ± 1.1) x 10"^, 
^ ; BiB+ cJ7r+) = {6.6tll ± 0-7) x 10~^ B{B r]K*^) = (19.8^1^ ± 1.7) x 10^^^ where the 
first error is statistical and the second systematic. For several other modes we report upper limits 
on their branching fractions; for example for the following flavor-changing neutral current decays, 
B{B Kl+l~) < 0.6 X 10-^ B{B K*l+l") < 2.5 x 10"^ at 90% Confldence Level (C.L.). 
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1 Introduction 



Flavor-changing neutral currents are forbidden at the tree level in the Standard Model (SM), hence 
such processes are possible only through penguin loops or suppressed tree amplitudes proportional 
to small couplings in hadronic flavor mixing (CKM matrix [l|]). These rare decays are interesting 
because their rates and kinematics are in principle sensitive to new heavy particles, predicted for 
example by supersymmetry models, which can enter the loop. Furthermore, it is possible to use 
some of these modes to search for direct CP-violation, measuring the CKM angle /? in a penguin 
environment and compare it with results at the tree level (from the charmonium modes). In 
the process, we can also test a great number of models just by measuring the branching fractions 
of penguin mediated B decays. 

The low rates of these decays and their large backgrounds make a high luminosity S-factory 
necessary for their searches. The results presented here are derived from 20.7 fb~^ delivered in 
1999 and 2000 by PEP-II g at the T{4S) peak and 2.6 fb"^ off peak (running at ~40 MeV 
below the resonance energy), for a total of ~22.7 million BB events. A detailed description of the 
BABAR detector can be found elsewhere Q. Charge conjugate states are assumed throughout and 
branching fractions are averaged accordingly. 

2 Analysis 

Much of the background in these rare decays can be reduced by exploiting the good charged particle 
identiflcation of BABAR. This is crucial in analyses involving B meson decays with no charm in 
the final state. dE/dx from the tracking devices provides good K/ti separation at p < 0.8 GeV/c, 
while the response of the internally refiecting ring imaging Cherenkov detector (DIRC) is excellent 
at higher momenta. Here only the measurement of the Cherenkov angle is used to discriminate K 
and vr (with better than 3 a separation up to 3.5 GeV/c), while for (lower momentum) composite 
particle daughters, so-called kaon selectors are used. These selectors combine information from all 
the relevant subsystems and have typical efficiencies of 80-90% with very low (~ 2-3%) vr — > i^T 
misidentification probability. The very tight selection criteria applied to electrons give an efficiency 
of ~88% with a corresponding vr — > e misidentification probability of ~0.15%. The typical muon 
selection efficiency for momenta greater than ~1 GeV/c is 60-70%, with vr ^ /x misidentification 
probability of ~2%. Performances of particle identification are tested on many control samples 
kinematically selected from data. Since the backgrounds are dominated by combinatorics from the 
continuum and since the continuum topology is jetty compared to the isotropic distribution of the 
signal, event shape variables are exploited to fight this kind of background. The main ones are a 
Fisher discriminant (already used at CLEO [^) and the cosine of the angle between the thrust axis 
of the B and the rest of the event. 

B meson candidates are selected using either the beam energy-substituted mass {mEs = 
\J {EexpY ~ {pbY, where is the momentum of the reconstructed B and E^xp the B candi- 
date expected energy) or the energy-constrained mass {mECi obtained via a kinematic fit of the 
measured candidate four momentum in the T(4S) frame with the constraint E"^ = -E'jeam)' 
the difference between the reconstructed B candidate and the beam energies (AE = E*^ — -f^^eam)' 
where the stars indicate variables evaluated in the T(4S) rest frame. The latter depends on the 
mass hypothesis, hence it is a good discriminant for different final states. All analyses use an un- 
binned maximum likelihood fit to extract the signal yields, while event counting methods are used 
just as cross checks. Most analyses select their cuts with their signal regions blinded. 
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Figure 1: Projections onto the variable niEs'- (a) — * (t^K^; (b) —> (j)K^\ (c) B'^ — > (f)K*~^] 
(d) B^ (f)K*^. In (c) the histogram is the sum of the two cjiK*^ channels, with the shaded area 
corresponding to the K*~^ — > K^tt^ decay chain. 



2.1 B cj)K^*^ 

In the B — > (pK^*^ modes we reconstruct the final states 4'K~^ , (pK^, (j)K*^ and (f)K*^ where 
the following intermediate states are recovered: K*^ K^tt^, K*^ —>■ K^t:^, K*^ K^it~ , 
(/) K^K^, K'^ TT^TT~ and vr^ 77. These decays are particularly interesting because 

they are dominated by 6 — > s{d)ss penguins, with gluonic and electroweak contributions, while 
other SM contributions are strongly suppressed p. In Figure || the mEs distributions for the 
selected candidates for all the previous modes are shown. The solid (dashed) line shows the PDF 
projection of the full fit (background only). The fits are mostly for illustration purposes, since the 
branching fraction results are derived, not from these, but from the maximum likelihood fits. A 
clear signal is visible in all channels and we report a first observation of the decays B^ (f)K*^ 



Decay Mode 


Branching Fraction 


S+ ^ (l)K+ 

B+ (pK*+ 
B^ 4)K*^ 
B+ (/)7r+ 


{7.7tiA ± 0-8) X 10-6 
(S.liy ±0.8) X 10-6 
(9.7l|^ ± 1.7) X 10-6 
(8.7l|f ± 1.1) X 10-6 
< 1.4 X 10-6@90% C.L. 



Table 1: Results of branching fraction measurements for the B (pK^*^ modes. 
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Figure 2: Projections of mEc for the B 
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mode. The projection is made by selecting events 



with signal likelihood (computed without mEc) exceeding a threshold that optimizes the expected 
sensitivity. 



Decay Mode 


Branching Fraction 


B+ ujK+ 

B+ UJTT+ 

B^ ^ ojK^ 

B^ ^ UJTT^ 


< 4 X 10'^ @90% C.L. 
{6.6tli ± 0.7) X 10-6 

< 14 X 10-6 @90% C.L. 

< 4 X 10-6 @90% C.L. 



Table 2: Results of branching fraction measurements for the B — > luX modes. 



and B^ —>■ cpK^. The final results for these modes are reported in Table ||. They are consistent 
with previously reported measurements [0] and, within errors, are consistent with isospin invariance 
under the assumption of penguin diagram dominance. 

2.2 B ujX 

B mesons are reconstructed also in their decays into ujn^ , ujt:^ , loK^, ojK'g, with the to decaying into 
three pions. A signal is found solely for the mode B — > uj-k^ (Figure and the branching fraction 
measurement reported in Table |2| is consistent with those previously reported |^]; in particular we 
find B{B^ — > W7r+) > B{B~^ — >• uK^)^ as expected. Improved upper limits are found for all the 
other decay channels. 

2.3 B r]K* 

We have also analyzed the modes B r]K*~^ and B — > t]K*^, where the t] decays into two photons, 
the K*^ into K^tt~, the K*~^ into Kg7r~^ and the Kg into tt^it~. The niEC projection plot is 
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Figure 3: B candidate invariant mass for rjK*^. The histogram represents data and the 

smooth curve the fit function. 



shown in Figure ^ for the r]K*'^ channel, where a clear signal is present and for which we measure 
a branching fraction of (19.8l5g ^ ^■'^) ^ 10^^. We also find an upper limit of 33.9 x 10^^ for the 
branching fraction of the r]K*~^ mode at 90% C.L. The results are consistent with those previously 
reported by the CLEO Collaboration Q and, at least in the r]K*^ mode, confirm their rather high 
branching fraction measurement. 

2.4 B K(*H+l- 

The dominant contributions for the processes B K(*h+l- in the SM are the so-called electroweak 
radiative penguins. These decays have been reconstructed in the modes where I = e, /i, the K*^ 
decays into K~^7r~ and the K*~^ into Kg{Kg — > 7r"'"7r~)7r"'". These processes have a very clean 
experimental signature, due to the presence of both a lepton pair and a kaon in the final state. 

The main challenge lies in understanding and characterizing the background. To avoid biases 
as much as possible, a blind analysis is performed, where not only the signal region, but also the 
sidebands used to determine the background and its normalization, are blinded. Several control 
samples from data are used to verify the reconstruction efficiencies as determined from Monte Carlo 
(MC) simulation. In particular B — > J/ipK^*^ and B tlj{2S)K^*^ events are used, since they have 
the same topology as the signal for these modes. Since these decays also constitute a dangerous 
background, quite complex cuts in the AE and the lepton pair mass plane are defined to exclude 
these events, taking into account bremsstrahlung effects; they are shown in Figure ^ where the 
shaded areas are vetoed. (The same veto is applied to the K* modes.) For reference, the two 
horizontal lines bound the region in which most signal events are found. The charmonium veto 
removes these backgrounds not only from the signal region, but also from the sideband region, 
simplifying the description of the background in the fits. 

No excess of events in the signal regions is observed, as can be seen in Figure |5|, where the AE, 
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Figure 4: The veto regions in the AE vs. mi+i- plane that are populated hy B ^ J/ij)K and 
B i;{2S)K events. 
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Figure 5: Scatter plots of AE vs uies after all analysis cuts in (a) B^ —i- K^e^e , (b) B^ —f 
K + ii+ii-, (c) B^ K*^e+e-, (d) 5° ^ K*^^+n-, (e) 5° ^ i^ge+e", (f) B^ K^sf^+fi', (g) 
5+ ^ i^^7r+)e+e-, (h) 5+ K*+{-^ K^7r+)n+n-. The small rectangles indicate the 

signal region, which is used only for optimizing event selection criteria. 
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Decay Mode 


Branching 
Fraction 
upper limits 
(@90% C.L.) 


B+ K+e+e- 
B^ K^ji+H- 


< 1.3 X lO"*^ 

< 0.9 X 10-6 

< 4.5 X 10-<^ 

< 4.7 X 10-6 


B Kl+l- 


< 0.6 X 10-6 


B^ iC*V+^- 
BO _^ K*^e+e- 
B+ K*+n+n- 
B+ ^ K*+e+e- 


< 3.6 X 10^6 

< 5.0 X 10-6 

< 17.5 X 10-6 

< 10.0 X 10^6 


B K*l+l- 


< 2.5 X 10-6 



Table 3: Results of branching fraction measurements for the B K^*^l'^l modes. 
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Figure 6: Comparison of AE shapes in several data control samples: (a) B^ J/ip{^ e^e )K~^ , 
(b) B+ J/V'(^ ii+ii-)K+, (c) B^ J/i^{^ e+e-)i^*o, (d) 5° ^ J/i^{^ ^+^-)K*^, (e) 
ijo ^ e+e-)_Kl (f) B^ ^ J/V(- (g) B+ ^ J/^(^ e+e-)K*+{^ KV), (h) 

B+ — > J/tp{-^ fj,'^fj,~)K*+{-^ Kg-ir~^). The normalization is absolute. 
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niES scatter plots are shown for all decay modes. Very low level of background is present. As a 
cross check, after removing the J/ip^ V'(2S) veto, we repeat the analysis on the J/ipll and ■ip{2S)ll 
control samples. Not only is a clear signal found, but the resulting AE shapes and normalizations 
are in very good agreement with the simulation, as shown in Figure |6|, where the points with error 
bars show the on-resonance data, and the solid histograms the predictions of the charmonium MC 
simulation. All the analysis selection criteria have been applied except for the charmonium veto, 
which is reversed. The results are shown in Table ^. The combined limits for B — > K^*^ll are now 
very close to the theoretical expectations for the SM [^] and represent a significant improvement 
over previous results n m pi. 



3 Summary 

With about 23 millions BB events BABAR has already very interesting results on penguin mediated 
decays. We have reported the first observation of the decays B~^ — > (f)K*^ and B^ — > (f)K^ for which 
we found the branching fractions (9.7;1;3 4 ± 1.7) x 10~^ and (S.ltfs ±0.8) x 10~^ respectively, and 
the upper limit set on the B{B — s- Kll), < 0.6 x 10~^, is very close to the theoretical expectations for 
the SM. Other branching fractions have been measured for penguin mediated B decays: B{B^ 
(j)K+) = {7.7tli ± 0.8) x 10-6, 5(^0 ^ ^K*'^) = (8.7t2-5 ± i.i) x 10"^, B{B+ ujtt+) = 
{G.Gtli ± 0.7) X 10-^ B{B rjK*^) = {W.StH ± 1.7) x 10^*^. 
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